Patterns on a Roll: A Method of Continuous Feed Nanoprinting 
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Exploiting elastic instability in thin films has proven a ro- 
bust method for creating complex patterns and structures 
across a wide range of lengthscales. Even the simplest of 
systems, an elastic membrane with a lattice of pores, under 
mechanical strain, generates complex patterns featuring 
long-range orientational order. When we promote this sys- 
^ tern to a curved surface, in particular, a cylindrical mem- 
^ brane, a novel set of features, patterns and broken sym- 
j metries appears. The newfound periodicity of the cylinder 
^ allows for a novel continuous method for nanoprinting. 
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Again and again, periodicity enables the design and descrip- 
tion of the materials we use to manipulate light, electrons, and 
other wave phenomena. Elastic instabilities have long been 
utilized as a means of generating patterns in thin films 1-11 . 
Generating, by self assembly, increasingly complex patterns 
on ever larger substrates with even finer structure remains a 
constantly receding goal. Mastery of such processes would 
revolutionize the fabrication and design of novel materials 
with specific properties. Recently, several groups have ex- 
ploited an elastic instability in a mechanically stressed elastic 
sheet perforated by a regular array of circular holes to gen- 
j erate intricate structures. Upon compression, the holes snap 
into elongated slits, revealing a diamond plate pattern fea- 
turing long-ranged orientational order 12-14 , a robust mech- 
anism acting from the macroscopic to the nanoscale. This 
long-range order can be faithfully transmitted to substrates as 
printed patterns of metallic nanonparticles, resulting in opti- 
cally interesting nanoscale materials. This method of one-step 
nano-assembly is limited by both the size of the elastic sheet 
and, more importantly, the limited number of possible sym- 
metries that can be cast: there are five Bravais lattices in two- 
dimensions. To bypass these limitations, here we report on 
compression induced patterns arising from circular holes on a 
cylinder. The additional symmetry afforded us by the period- 
icity of the cylinder combined with the long-range elastic in- 
teractions between the slits leads to an extended library of mo- 
tifs. By capitalizing on the periodicity of the cylinder, one can 
upgrade from the "sheet-at-a-time" contact printing method 
for flat membranes to continuous-feed printing. Thereby pre- 
viously unachievable patterns can be transferred to arbitrarily 
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long flat substrates with possibly interesting optical properties. 

The collapse of holes in an elastic membrane - both flat and 
curved - occurs in the highly non-linear regime of elasticity, 
complicating the analysis. Technological advances have made 
computational techniques, for instance finite element simula- 
tions, increasingly feasible. Only finite element simulations 
using specific models of nonlinear elasticity capture the entire 
process of the holes collapsing 14,15 . These nonlinear effects 
are needed to predict the final shape of the collapsed holes 
and not merely their orientation. As the complexity of the sys- 
tem grows, such methods tend to be more accurate but lack 
the ability to distinguish between mechanisms responsible for 
different phenomena. Analytic calculations generate intuitive 
solutions and can form the basis for subsequent numerical cal- 
culations. With the sole assumption that each hole collapses 
to some elongated shape, our model uses only linear elasticity 
to successfully predict the orientational order in the diamond 
plate pattern and the herringbone pattern formed from an un- 
derlying triangular lattice. Not only does our model shed light 
on the interactions in the system, it greatly facilitates the ra- 
tional design of other patterns and devices. 

Taking the theory of cracks as the starting point, the far field 
deformation of a collapsed hole can be approximated by a con- 
tinuous distribution of parallel dislocations. The details of the 
shape of the hole, and thus the exact distribution of disloca- 
tions, result from nonlinearities in the elasticity. Since we con- 
sider only the final elongated slit, the first Fourier mode, a dis- 
location dipole of strength b with dipole vector d described by 
the Burgers vector b(r) = z x db ( - 8(r - d/2) + 8(r + d/2)) , 
sufficiently captures the deformation of a single hole. Using 
linear elasticity theory, the interaction energy between a pair 
of dipoles with dipole vectors di and d2 centered at r = and 
r = R, respectively, is 



E = — 



Y 2 b 2 did 2 
TlR 2 



cos(0i + e 2 ) sin(6i) sin(0 2 ) + 



(1) 



where Y 2 is the two dimensional Young's modulus and 
cos(0i) = R di and cos(02) = R d 2 are the angles each 
dipole makes with the vector connecting them. For any ar- 
rangement of holes, we simply minimize the total energy - 
the sum of all pairwise terms for a given set of holes - for 
each of the dipole angles to find the resulting groundstate con- 
figuration upon hydrostatic compression. 

The introduction of curvature immediately complicates the 
system. Not only must the elastic energy for the membrane 
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Fig. 1 (A) A cylinder of radius R with two dislocation dipoles di and d2. (B) The unrolled version of the cylinder with the first two image 
cylinders at 2nnR% n = ±1. 



change to include curvature terms, the concept of a lattice is 
no longer well defined. By adopting the topology of a cylinder, 
we may bypass these complications. Because the cylinder is 
isometric to the plane, it has zero Gaussian curvature. We 
may neglect additional bending energy if the thickness of the 
membrane is much less than all other relevant length scales. 

In order to calculate the interaction between two holes on 
a mechanically compressed cylindrical membrane, we begin 
with the pairwise interaction energy, Eq. (1), of two holes 
on a flat membrane and consider the effect of the additional 
periodicity given by the cylinder. 

Our geometry consists of an infinite elastic cylindrical shell 
of radius R with its axis parallel to z passing through the ori- 
gin with two dislocation dipoles on it. Imagine cutting the 
cylinder along the line r = 7?, = 0, and unrolling it onto the 
plane such that x = 7?0, y = z. The two dipoles are now lo- 
cated at di = {x\ , y\ } and d2 = {^2,^2}- They interact not only 
through the shortest path and the first two replicas located at 
di ± 2nRx, as they would if one simply considered periodic 
boundary conditions, but, due to the range of the potential, the 
interaction must also include terms from each of the infinite 
number of "image" holes. Note that the calculation on a cylin- 
der is the same as the calculation on an infinite elastic sheet 
with a copy of the unrolled cylinder located every 2nnR along 
the x direction (see FIG. 1). Therefore, the reduced interac- 
tion energy, *£ = ER 2 / faith 2 d\d,2) , between two dislocation 
dipoles on a cylinder is given by, 



be summed, 



cos(0i + e 2 - 26 n ) sin(6i - Q n ) sin(6i - 9„) 



n^oo (X + 2nn) + Y 2 

(2) 

where X = (x 2 -x 1 )/R,Y = (y 2 -yi)/R, and tan0„ = Y/(X- 
2izn) is the angle between dipole 1 and the n th image of dipole 
2. Upon expanding the trigonometric functions in equation 
(2), the reduced energy consists of 5 terms, each of which may 



00 4 Y^~^(v 

1 1X61,62) 1 



-2nn) 



[(X + 2nn) 2 + Y 2 ] : 



(3) 



where the fl* (61,62) are given in Table 1. We perform the 
tedious but straightforward sum in the appendix and derive 
the final expression for the energy. 




Fig. 2 The simple model using linear elasticity theory in equation 
(3) faithfully reproduces observed results for the square lattice (data 
on left). 

The effect of long ranged elastic interactions serves to rein- 
force the same angular dipole interaction seen in the flat case 
for dipole pairs located at the same height on the cylinder, 

e(x,y = o) = - Y 2^d 2 ^ cos (e 1 + e 2 ) S i n e lS i n e 2< 

Wrapping a lattice around a cylinder breaks the translational 
symmetry of the plane along one direction. The orientation 
and magnitude of this periodic direction with respect to the 
lattice vectors and spacing, give new degrees of freedom with 
which to control patterns. To gain a handle on this new phase 
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Table 1 Expressions for the functions %(0i , 62) used in the sum (3). 
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cos(0i + ©2) sin0i sin 02 



space of possible patterns, we examine only those resulting 
from achiral lattice wrappings. 

Before examining the plethora of complex lattices, it be- 
hooves us to consider the simple example of a square lattice, 
the results of which have been verified by a physical model 
(see Fig. 2). On a flat membrane, it is well known that a square 
lattice will produce a diamond plate pattern. Although the pat- 
tern produced by the equivalent square lattice on a cylinder 
resembles a diamond plate pattern, the orientations of the col- 
lapsed holes are neither orthogonal to one another, nor do they 
lie along principal lattice directions. 

Fig. 3 demonstrates the effect of changing the wrapping di- 
rection of simple lattices. The triangular lattice (Fig. 3 A) pro- 
duces stripes of alternating angles when one of its lattice vec- 
tors is aligned with the wrapping direction, and it produces a 
compressed diamond plate pattern when parallel to the cylin- 
drical axis. The honeycomb lattice, when wrapped in the two 
aforementioned directions (Fig. 3B), produces very different 
results. Of particular interest, the latter wrapping breaks chiral 
symmetry when the lattice is an odd number of hexagons tall. 

A systematic study of variations of these simple lattices 
yields very precise control over the exact orientations of the 
dislocation dipoles for a few classes of patterns. As variations 
on a theme, they can only add incrementally to our library of 
patterns. We must turn to increasingly complex lattices in or- 
der to broaden the range of accessible patterns. 

The kagome lattice represents a perfect case study, as it has 
recently gained popularity in the condensed mater literature 
for its novel mechanical and vibrational properties 16-18 . The 
kagome lattice consists of a triangular lattice ai = {2a, 0} a 2 = 
{a,V3a} with a basis of three points bi = {0,0}, b 2 = 
{<2,0}, b3 = {|, y -^ L }. The resulting pattern retains remnants 
of the symmetry of the underlying lattice, as shown in Fig. 4. 

Varying the relative alignments of the lattice and cylindrical 
axis breaks the lattice symmetry and introduces a new degree 
of freedom to generate patterns. Moreover, unlike the dia- 
mond plate lattice on a flat membrane, the patterns on cylin- 
ders are sensitive to the number of holes both in circumference 
and in height. This immediately raises the question: how do 
we generate extensive patterns? Because cylinders are isomet- 
ric to the plane, Gaussian curvature does not hinder the ability 
to transfer patterns onto a flat substrate. By capitalizing on the 




Fig. 3 The resulting patterns depend not only on the type of lattice 
but on its alignment with respect to the cylindrical axis, for example 
numerical results for triangular (A) and honeycomb (B) lattices. 



periodicity of the cylinder, one can upgrade from the "sheet- 
at-a-time" contact printing method outlined in 13 to continu- 
ous transfer. By filling the interior of the cylinder with an ink 
containing, for instance, nanoparticles, quantum dots, or poly- 
mers, the exact orientation and position of the holes may be 
transferred to a flat substrate of any dimension, shown in Fig. 
4. The newly patterned surface adopts the same symmetries 
and optical properties of the cylindrical lattice not previously 
achievable with a flat membrane. Likewise, by varying the ink 
pressure, changes in the pattern spacing and slit geometry can 
be controlled on the fly. Once again, periodicity makes this 
possible. 
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Fig. 4 An example application uses the pattern generated by the 
square lattice on a cylinder as a rolling printer onto a flat substrate. 



A Calculation of Geometric Sums 

In order to calculate the five sums in equation (3), we first 
~ 1 sinhg 



note that j^ + 2)(n) 2 + 9 2 

The five sums are given by: 
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where H = cosX — cosh Y. 
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Combining these we arrive at a final, albeit complicated, ex- 
pression for the pairwise interaction of two dislocation dipole 
on a cylinder, 

Y 2 b 2 did 2 n r2# 2 sinh7 
E ~ 4^fP [ Y 

+ 2if sinXsinh7(sin2ei +sin2e 2 -sin2(ei +e 2 )) 

+ (3H + cos X cosh 27 -cos 2X cosh Y) 

x (cos20i +cos20 2 - cos 2(0i +62)) 

+ Y ((cos2X - 3) sinhF + cosX sinh27) cos2(0i + 6 2 ) 

+ Y ((cosh27 - 3) sinX + coshF sin2X) cos2(0i + 9 2 ) 
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